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Abstract 

We study the inclusion of running coupling corrections into the non-linear small-x 
JIMWLK and BK evolution equations by resumming all powers of agNf in the evolu- 
tion kernels. We demonstrate that the running coupling corrections are included in the 
JIMWLK/BK evolution kernel by replacing the fixed coupling constant in it with 

"'^^^afijit'^)^^^ , where ri and r2 are transverse distances between the emitted gluon and 
the harder gluon (or quark) off of which it was emitted to the left and to the right of 
the interaction with the target. In the formalism of Mueller's dipole model ri and r2 are 
the transverse sizes of "daughter" dipoles produced in one step of the dipole evolution. 
The scale i? is a function of two-dimensional vectors ri and r2, the exact form of which is 
scheme-dependent. We propose using a particular scheme which gives us R as an explicit 
function of ri and r2. 



1 Introduction 



In the recent years there has been a lot of progress in small-x physics due to developments in 
the area of parton saturation and Color Glass Condensate (CGC) [1-21]. Among other things 
the CGC led to a new way of calculating the hadronic and nuclear structure functions and total 
cross sections in deep inelastic scattering (DIS) at small values of Bjorken x variable. Accord- 
ing to the CGC approach to high energy processes, one first has to calculate an observable in 
question in the quasi-classical limit of the McLerran-Venugopalan model [3-5] which resums 
all multiple rescatterings in the target hadron or nucleus. After that one has to include the 
quantum evolution corrections resumming all powers of ag In I/xbj along with all the multiple 
rescatterings. Such corrections are included in the general case of a large target by the Jalilian- 
Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) functional integro-differential 
equation [9-16], or, if the large- A'^c limit is imposed, by the Balitsky-Kovchegov (BK) integro- 
differential evolution equation [17-21] based on Mueller's dipole model [22-25]. The JIMWLK 
and BK evolution equations unitarize the Balitsky-Fadin-Kuraev-Lipatov (BFKL) linear evo- 
lution equation [26,27]. For detailed reviews of the physics of the Color Glass Condensate we 
refer the reader to [28-30]. 

Both the JIMWLK and BK evolution equations resum leading logarithmic Og Inl/xsj cor- 
rections with Og the coupling constant. At this leading order the running coupling corrections 
to the JIMWLK and BK evolution kernels are negligible next-to-leading order (NLO) correc- 
tions. A running coupling correction would bring in powers of, for instance, Inl/xsj, which 
are not leading logarithms anymore. Hence both JIMWLK and BK evolution equations do not 
include any running coupling corrections in their kernels. The drawback of this lack of run- 
ning coupling corrections is that the scale of the coupling constant to be used in solving these 
evolution equations is not known. Indeed, as was argued originally by McLerran and Venu- 
gopalan [3-5] and confirmed by the numerical solutions of JIMWLK and BK equations [31-34], 
the high parton density in the small-z hadronic and nuclear wave functions gives rise to a hard 
momentum scale — the saturation scale Qs- For small enough x and for large enough nuclei 
this scale becomes much larger than the QCD confinement scale, Qs ^ Aqcd- The existence 
of a large intrinsic momentum scale leads to the expectation that this scale would enter in the 
argument of the running coupling constant making it small and allowing for a perturbative 
description of the relevant physical processes. However, until now this expectation has never 
been confirmed by explicit calculations. 

In the past there have been several good guesses of the scale of the running coupling in 
the JIMWLK and BK kernels in the literature [33,35]. A resummation of all-order running 
coupling corrections for the linear BFKL equation in momentum space was first performed by 
Levin in [36] by imposing the conformal bootstrap condition. There it was first observed that 
to set the scale of the running coupling constant in the BFKL kernel one has to replace a single 
factor of as by the "triumvirate" of couplings asag/c^s with each coupling having a different 
argument [36]. 

In this paper we calculate the scale of the running coupling in the JIMWLK and BK evolu- 
tion kernels. Our strategy is similar to [37]: we note however, that [37] relies on the dispersive 
method to determine the running coupling corrections, while below we use a purely diagram- 
matic approach. We concentrate on corrections due to fermion (quark) bubble diagrams, which 
bring in factors of UgNf. Indeed some factors of Nf may come from the QCD beta-function 
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(see Eq. (j^)) below), while other factors of Nf may come in from conformal (non-running cou- 
pling) NLO (and higher order) corrections [38-40]. While we do not know how to separate 
the two contributions uniquely, we propose a way of distinguishing them guided by UV diver- 
gences. This leaves us with an uncertainty with respect to finite contributions in separating the 
conformal and the running coupling factors of Og Nf that influence the scale of the obtained 
running coupling constant in a way reminiscent of the scheme dependence. Once we pick a 
certain way of singling out the factors of as Nf coming from the QCD beta-function, we replace 
Nf ^ —6% j32 { "completing" Nf to the full beta-function) and obtain all the running coupling 
corrections to the JIMWLK and BK kernels at the one-loop beta-function level. 

The paper is structured as follows. We begin in Section |21 by calculating the lowest order 
fermion bubble correction to the JIMWLK and BK kernels, as shown in Figs. and |2l in the 
framework of the light cone perturbation theory (LCPT) [41,42]. We note that the diagrams 
in Figs. andEK' give a new kind of evolution kernel, which does not look like a higher order 
correction to the LO JIMWLK or BK kernels. We analyze the problem in Section IHl where we 
propose a subtraction procedure to single out the part of these diagrams' contribution giving 
the running coupling correction. There we show that this subtraction procedure is not unique 
and introduces a scheme dependence into the scale of the running coupling. 

In Section|3]we resum fermion bubble corrections to all orders, and, after the Nf —>■ —6 vr /?2 
replacement obtain the JIMWLK evolution kernel with the running coupling correction given 
by Eq. ()88|) in transverse momentum space as a double Fourier transform. The corresponding 
BK kernel is obtained from Eq. ()88|) using Eq. ()90|). Notice that the running coupling comes in 
as a "triumvirate" originally derived by Levin for the BFKL evolution equation [36]. Fourier- 
transforming the running couplings into transverse coordinate space is more involved since one 
encounters integration over Landau pole leading to power corrections. A careful treatment of 
the uncertainties associated with power corrections in small-x evolution was performed in [37]. 
Here we calculate the Fourier transforms by simply ignoring those corrections and by using the 
Brodsky-Lepage-Mackenzie (BLM) method [43] to set the scale of the running coupling. The 
JIMWLK kernel with the running coupling corrections in the transverse coordinate space is 
given in Eq. (PH|) . 

We conclude in Section El by explicitly writing down the full JIMWLK Hamiltonian with the 
running coupling corrections in Eq. ()99|) and the full BK evolution equation with the running 
coupling corrections in Eq. (llOlj) and by discussing various limits of the obtained result. 

We note that our analysis is complimentary to [37], where the running coupling correction to 
the JIMWLK and BK kernels was determined using the dispersive method. Our result for the 
all-order series of A^j-terms is the same as in [37]. However, using the diagrammatic approach, 
we have been able to identify the structure of that series as coming from a "triumvirate" of the 
coupling constants in Eq. ()98|) . which is an exact result in the transverse momentum space and 
a better approximation of the full answer in the transverse coordinate space. 

2 Leading Order Fermion Bubbles 

Our goal in this work is to resum all OgNf corrections to the leading logarithmic non-linear 
JIMWLK and BK small- a; evolution equations [8-21] (for review see [28-30,44]). After extract- 
ing the running coupling a^A^j-corrections out of all possible asNf terms, the complete running 
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coupling correction to the JIMWLK and BK evolution kernels would then be easy to obtain by 
replacing 

Nj-^-Qnf52 (1) 

in the former, where 




ABC 



Figure 1: Diagrams giving the leading UgNj correction to the kernels of JIMWLK and BK 
small-x evolution equations. The thick dots on gluon and quark lines denote interactions with 
the target. 




A' C 



Figure 2: Diagrams with instantaneous parts of gluon propagators giving the leading UgNj 
correction to the kernels of JIMWLK and BK small-x evolution equations. There is no analog 
of Fig. Ql B. All the lines are implied to he labeled in the same way as in Fig. 

To resum UsNf corrections we begin by considering the lowest order diagrams for one step 
of small-x evolution containing a single quark bubble. These diagrams give the lowest order 
asNf correction to the JIMWLK and/or BK evolution kernels and are shown in Fig. ^ The 
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diagrams are time-ordered as they are drawn according to the rules of LCPT [41,42]. Gluon 
hnes in Figs. and C also have instantaneous/longitudinal counterparts [41,42], shown 
in diagrams A' and C in Fig. |21 (The virtual gluon on the left side of Fig. ^3 can not be 
instantaneous, since the produced gluon on the right of Fig. CJ3 can only be transverse and a 
longitudinal gluon can not interfere with a transverse gluon, as will be seen in the calculations 
done below.) 



2.1 Diagrams A and A' 

To calculate the forward scattering amplitude in Figs. ^A. and|2lA.' we first need to calculate the 
wave function of a dipole emitting a gluon which then, in turn, splits into a quark-anti-quark 
pair, i.e., the part of the diagrams A and A' located on one side of the interaction with the 
target. The calculation is similar to what is presented in [45]. We will work in the = 
light cone gauge in the framework of the light cone perturbation theory [41,42,46,47]. The 
momentum space wave function of a dipole (or a single (anti-) quark) splitting into a gluon 
which in turn splits into a qq pair with the transverse momenta ki and ^2 of the quark and the 
anti-quark with the quark carrying a fraction a of the gluon's longitudinal ("plus") momentum 
is [45] 



¥Jl^^{k^,k,,a) = [r]e„.® [t1fV^i^,l.,(fci,fc2,a) = -2/[r]em® [tif 



A=±l 



2 



^2 



e*^ • (fci + fca) 
{ki + k 

^ kl{l-a) + kla [t]en.®[t]f^2(i_^) + fc2^- (3) 

Here A = ±1 is the internal gluon's polarization: the gluon polarization vector for transverse 
gluons is given by = (0,0, e'^) with e'^ = (1 + iX)/\^. The instantaneous diagram from 
Fig.|2K' gives the second term on the right hand side of Eq. The produced quark and anti- 
quark are massless, which is sufficient for our purposes of determining the scale of the running 
coupling. (Ti = ±1 and (T2 = ±1 are quark and anti-quark helicities correspondingly (defined 
as in [45]). The fraction the of gluon's "plus" momentum carried by the quark is denoted by 
a = ki+/{ki+ + k2+). The wave function also contains a color factor [t"]em ® [t"]f consisting of 
two color matrices originating in the quark-gluon vertices at the points of emission of the gluon 
and its splitting into a qq pair. 

It is convenient to rewrite Eq. Q in terms of a different set of transverse momenta. Defining 
the momentum of the gluon q = ki + k2 and k = ki{l — a) — k2a, and noting that fc^(l — a) -|- 
^20! = fc^ + q^a(l — a), we write 

e*'^ ■ q e'^ ■ fc (1 — 2a + Aai) (5o-io-2 ^ 2 a {1 — a) 60 

x=±i 

Performing the summation over gluon polarizations A yields 



^' ^ + q^a{l — a) k^ + q^a(l — a) 



(5) 
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where denotes the ith component of vector q and the sum over repeated indices i, j = 1, 2 is 
imphed. Here ei2 = 1 = —£21, Cii = £22 = 0, and, assuming summation over repeating indices, 

To find the contribution of the diagrams in Figs. ^A. and to the next-to-leading order 
(NLO) evolution kernel we first have to transform the wave function from Eq. (0) into transverse 
coordinate space 



(1) 



'^^^1 ^^^2 ^-ifci-(zi-cc„)-ifc2-(z2-3:„) ^(1) 



(2 7r)2 (2 7r) 



[fci,fc2,a). (6) 



Here the transverse coordinates of the quark and the anti-quark are taken to be Zi and Z2 
correspondingly. The gluon in Fig. ^ can be emitted either off the quark or off the anti-quark 
in the incoming "parent" dipole. The transverse coordinates of the quark and the anti-quark 
in the "parent" dipole are Xq and Xi. In Eq. ^ we labeled them Xm with m = 0, 1 depending 
on whether the gluon was emitted off the quark or off the anti-quark. 
In terms of transverse momenta k and q Eq. (jH)) can be written as 



(1) 

0-1,0-2 



Xyyi^ Zi Xjji, <yj 



where 



and 



(2 7r)2 {2 Try 



Z12 — ^1 — ^2 



(7) 



a Zi 



[1 - a) Z2 



(9) 



is the transverse position of the gluon. 

Substituting the wave function from Eq. (0) into Eq. ((7j) and performing the integrations 
over k and q yields 



1 



(Z - Xr. 



{2n) 



a (1 — a) z 



12 



X [(1 - 2a)6ij + icTieij] 6, 



^12 



010-2 ^2 



^9' 



a{l — a) 6a 



12 



(27r)2 (z — Xm)"^ + a{l — a) z 



2 ■ 

12 



(10) 



To calculate the diagram in Figs. ^A. and |2K' using the wave function from Eq. (fTUI) in 
a general case we have to include the interaction with the target by defining path-ordered 
exponential factors in the fundamental representation 



P exp 



-ig 



0,a^) 



fill 
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With the help of Eq. (jllj) we can write down 





+ 



d'zid^Z2 al /CfL°(a^o, xi; z,, z^) U^^t^ ® Ul/ 2 tr(t^f/,,t"f/tj ln(l/xBj), 

(12) 

where xsj is the Bjorken x variable. In arriving at Eq. p2|) we have defined the NLO contri- 
bution to the JIMWLK kernel coming from the diagrams in Figs. ^A. andEK', labeled /Cf^°, 
by multiplying the wave function in Eq. (fTUI) by its complex conjugate, summing the obtained 
expression over the helicities of the quark and the anti-quark in the produced pair and over Nf 
quark flavors, and integrating over a: 



«M ^^^^{Xo,Xi] Zi,Z2) 



Nf 



2(4 7r) 



da 



E 







0'1,0"2 = 



X V, 



(1) 



^Zi - Xq, Z2 



[zi - Xi,Z2 - xi,a} 



(13) 



In this definition of /Cf we use the wave function tjj^-^^ 0-2 , which is different from the full wave 

by the fact that the color matrices are included in ^^^^,0-2 and are 



function ^a^aa from Eq. 
not included in ^/'ailo-j- We have used 'ip^i,a2 to define the JIMWLK kernel /Cf'"^ because the 
color matrices were already included in the forward amplitude in Eq. p2|) . A factor of 1/2 
was inserted in Eq. (jl3p to account for the factor of 2 introduced in the definition of /Cf"^ in 
Eq. (Ha). 

Substituting ip^^^ from Eq. (jlUj) into Eq. (jl3j) and summing over quark helicities yields 



NLO, 



Xo,Xi] Zi,Z2 



da 



X 



'1 - 2afzi2 ■ {z - xo) Z12 



xi) + eij{z - Xo)i{zi2)j eki{z - Xi)k{zi2)i 



.2 \2 
'12> 



-2a{l-a) (l-2a) 



[{z - XqY + a{l-a) ZI2] [{z - xiY 

Z12 ■ (Z - Xq) + Z12 



a(l — a) 



12J 



[Z-Xi, 



-12 



[{z - XqY + a{l-a) ZI2] [{z - xi] 



ail — a) 



12J 



+ 



4^2(1 -a) 



[(z — aso)^ + a (1 — a) z\.^ \{z — x^'^ + « (1 — «) 



(14) 



The integral over longitudinal momentum fraction a, while straightforward to perform, would 
not make the above expression any more transparent. When squaring from Eq. (HUD one 
gets a cross-product between the first and the second terms on the right hand side of Eq. (fTTH). 
given by the second term in the square brackets of Eq. ()14p . Terms like that are also present 
in other physical quantities, such as the qq production cross section calculated in [45]. 
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To obtain the contribution of the diagrams in Figs. ^A. and to the BK evolution kernel 
Kf^'^ we have to sum the wave function ^l/i-^.o-a over all possible emissions of the gluon off the 
quark and off the anti-quark, multiply the result by its complex conjugate, sum over quark and 
anti-quark helicities and Nf quark flavors, take a trace over color indices averaging over Nc 
colors of the incoming dipole and integrate over a 



N 



2 TV-NLO/^ nr. ■ ^ ^ \ ^"Z 



r ' ' 

n ai,a2=~i- m,n=0 



X 



[^alaii^l -Xo,Z2- £Co, a) ^1^2*^,(21 - Xi,Z2- Xi,a 



(15) 



(Note that the capital K denotes the kernel of the BK evolution equation, while the calligraphic 
/C is reserved for the JIMWLK evolution kernel.) Using the first line of Eq. (jSj) along with 
Eq. in Eq. (fTK|l one can show that 



m,n=0 



(16) 



For the reasons which will become apparent momentarily, it is more convenient to leave 
/Cf^*^ written in terms of integrals in transverse momentum space. Using Eqs. ©, (|7j) in 
Eq. ()1H|1 and summing over quark helicities yields 

1 

[xo, x,,z,, Z2) 4 iV; J aa J ^^^^^ ^^^^^ ^^^^^ ^^7^)2 ^ 



1 (1 - 2aYq ■ k k' ■ q' + q ■ q' k - k' - q - k' k - q' 



q2qi2 



+ 



+ 



k^ + q'^a{l — a) k'"^ + q''^a{\ — a] 
2a(l-a)(l-2a) [k q , k' ■ q' 



k^ + q^a{l — a) fc'^ + q''^a{l — a) 



j'2 



k^ + q'^a{l — a) k'"^ + q''^a{l — a 



(17) 



where we have used the identity 

eij qi kj ekiq'kk'i = q ■ q' k ■ k' - q ■ k' k ■ q' . 



2.2 Diagram B 

Unlike the diagram A, the diagram B in Fig. ^ looks more like a "typical" running coupling 
correction to the leading order JIMWLK/BK kernels. The contribution of the diagram B along 
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with its mirror-reflection with respect to the hne denoting the interaction with the target can 
be written as 




d'z «J/Cf O(cco,a;i;2) U^X^Ui/ Uf ln(l/xBj) (19) 



with the corresponding NLO contribution /C^^*^ to the JIMWLK kernel calculated using the 
rules of the light-cone perturbation theory [41,42]. We first decompose the kernel /C^'"^ into a 
sum of the contributions of the diagram Fig. HP (denoted /C^i^°) and its mirror-image (denoted 



t-NLO V 
'^2 right/- 



-NLO 



NLO 



(20) 



Below we will only calculate JC2ieA- construct /C^^^^^ one only has to replace a^o ^ in its 
argument. A simple calculation along the same lines as the calculation of the diagram A done 
above yields 



o^l ^2Mt{xo,Xi;z) 



Nf 



(4 7r 



da 



a{l-a) J (27r)2 (27r)2 (27r)^ 



e*^-qe^-k{l-2a + Xai)6^,^^ e*^' ■k{l-2a + AVi) e^' • q' 

(Ti,(T2 = -l A,A' = -1 ^ 



q'^a{l — a) 



with all the notation being the same as in the case of the diagram and A' = ±1 the 
polarization of the gluon interacting with the target. Different from /Cf^*-*, the kernel JC^Mt i^ 
Eq. ()21|) has a part of the color factor included in it: it includes of 1/2 coming from the color 
trace of the quark loop, which is required by the definition of /C^^*"* in Eq. (jl9j) . Similar to the 
above, to obtain the corresponding correction to the BK evolution kernel, we use 



(22) 



m.n=0 



[Similar relationships holds for for both /C^^Jft ^2'right separately.] 

We can simplify Eq. First we sum over the quark helicities and gluon polarizations to 
obtain 



^NLO 



2 left 1^0; ^li ^) 



da 



df^^d^ -.q.(.-.o)W{.-l) 



a{l-a) J (27r)^ (27r)2 (27r)^ 



X 



;i - 2ayq kk q' + q q' -q kk q' 



{q^fq 



12 



fc2 



q'^a{l — a) 



(23) 



The integral over k is UV- diver gent, as expected. We will regularize it by using dimensional 
regularization, for which purpose we have replaced d'^k/{2'KY d'^k/{27rY in Eq. ()22|) with d 
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the number of dimensions. Anticipating the integration over the angles of the vector k we also 
replace 



fc2 

kikj —5ij (24) 



in Eq. (j2Sl), obtaining 



1 

a{l-a) J (27r)2(27r)2 [q^f q'^ 





f T^T-o ^-7 \(l-2a? + d-l\ . (25) 



Now the A;-integral is easily doable (see, e.g., [48]) yielding 



1 





X 



2 (47r)'^/2 

Writing d = 2 — e and expanding around e = we get 

1 



a{l-a) J (27r)2(27r)2 (q^)^ q'^ 

r(-0 [qMl-c^)f" [{l-2af + d-l]. (26) 



l^,,,,,[x,,x,,z) - J da J e 



X |[(l-2a)2 + l] 



In ^2 + 7 - In 47r - 1 



where we replaced 1/e with Iu/xms- Integrating over a we obtain 

^'•^^^^ °' ^' ^ " Stt 7 (2vr)2(27r)2" In ^ 



(27) 



(2^ 



with = /i^g 47r e ^. 

While Eq. ()28j) is sufficiently simple for our later purposes, we can further simplify it by 
Fourier-transforming it into transverse coordinate space. A straightforward integration yields 
the NLO contribution to the JIMWLK kernel coming from the diagram B 



)C^I;^{x„ X,; z) = ^^ . In i-^ - I - 27 ^ (29) 

67r-^|2-a;or \z - x^Y y \z - Xq^ 3 

Similarly one can show that 



^St(^o, xr, z) = f ^ ■ ^ <; In — ^ - i - 27 !> (30) 



Nf z — xq z — xi \ 4 5 
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and 



«(xo, X,; z) = JC^tfti^o, X,; z) + /C^.^g^ t(^o, x,- z) 



Nf z-xq z-xi I, 4e"5-2'^ , 4e-5-2T 

111 -; 7^ h m ■; 



6 vr^ \z — a^oP \z — Xi\ 



\z — ajop IJ? 



MS 



z — a^ip /i^ 



(31) 



MS 



The corresponding contribution to the NLO BK kernel can be easily obtained from Eq. 1)311] 
using Eq. 

Recalling that the leading order (LO) JIMWLK kernel is given by 



}C^'^{xo,Xi; z) 



1 Z — Xq Z — Xi 



TT^ \Z — £CoP \z — iClP 



(32) 



we immediately see that adding /C^^*-* from Eq. (p?T|l to it yields 
}C^^{xo, Xi; z) + al /C^^°(a;o, Xi] z) 



m'^S l^^O'^l' - ^2 _ ™ 12 ■ U _ ^ 12 



TT^ \Z — Xq\^ \Z — Xi\ 



6 TT 



In 



4e-f-2^ 



+ ln 



4e-f-2^ 



(33) 



Anticipating the appearance of the full QCD beta-function we perform the replacement of 
Eq. dH) in Eq. ^ to obtain 



a,/C^«(^o,^i;^) + «J/Cf°(a;o,^i;^) = ^ " 



X <( 1 - /32 



In 



4e-f-27 4e-f-2^ 
1- In 

19 2 ' I 19 2 

^-^Ol /^MS I^-=^iIVmS 



(34) 



Now one can readily see that the diagram B in Fig. ^ gives a contribution to the one-loop 
running coupling correction to the LO JIMWLK and BK kernels, as expected. 

2.3 Diagrams C and C 

The contribution of the diagrams in Figs. and ^P' along with their mirror-reflections can 
be written as 




+ 




j d^z al}C^'^''{xo,x,;z) U^^ ® U^t^ MVa^Bj 



(35) 



where z is still the gluon's transverse coordinate which we choose to keep explicitly even though, 
since both gluon lines are now completely virtual, they do not interact with the target. Instead 
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of calculating the NLO correction to the JIMWLK kernel AT^^*^ coming from the diagrams C 
and C explicitly we will use the conservation of probability condition, which states that, in 
the absence of interaction with the target, the sum of all three diagrams in Fig. ^ (along with 
the mirror images of the diagrams B and C reflected with respect to the line representing the 
interaction with the target) gives zero. Intuitively this condition is clear: in the absence of 
interactions there will be no contribution to the evolution kernel. In the diagrammatic sense, 
adding up all the graphs in Fig. ^ corresponds to summing over all the cuts for the diagram 
of the gluon emission with a quark bubble correction. Similarly, if the interactions are absent, 
the sum of the diagrams in Fig. |21 along with the mirror- reflection of the diagram C also gives 
zero. This probability conservation condition was originally used by Mueller to calculate the 
virtual correction to the leading order gluon emission in the dipole evolution kernel in [22]. In 
our case it formally reads 

j dh^(fz2]C^^^{xo,Xi;Zi,Z2) + J IC^^^ixo,xr,z) + J dh }C^^^{xo,Xi; z) = 0. (36) 



Noting that 



d'^Zi d^Z2 = d^z d^Zi2 



(37) 



with Zi2 and z defined in Eqs. (jS} and Q above, an explicit diagram calculation (keeping 
all the transverse momenta fixed in momentum space) yields an even stronger identity than 

Eq. (jnni): 



K^'^'^ixo^x.-z) 



j d2^i2/CfLO(a;o,a;i;zi,Z2)- /Cf ^(xo, z). 



(3^ 



^NLO is given explicitly in Eq. (jHH). Using the momentum-space expression (fTTj) 

for X:f we get 



1 (1 - 2ayq kk q' + q q' -qkk q' 



q2qi2 



+ 



k^ + q'^a{l — a) fc^ + q''^a{l — a) 



k^ + q^a[l-a) 



fc2 + q'^a{l - a) 



(39) 



where the fc-integral is UV- divergent, which we regularize using dimensional regularization. 
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With the help of Eq. we rewrite Eq. as 

1 



^•^'^ ^ [(l-2a)^ + ci-l] 



+ 4«2^1-«)2 j 



1 



k'^ + g2a(l — a) k'^ + q''^a{l — a) 



(40) 



where we put c? = 2 in the second term in the curly brackets since the integral in that term is 
not divergent. Performing the fc-integrals yields 

1 

r I- A'i, 

j2 ^ i^NLO / 



d ZulCi {xo,Xi; Zi,Z2) = ANf I da 



(27r)2 (27r)2 



qq'l 1 



q^q'^d (4 7r)'^/2 



r 1 



[a{l-a)f^-^ [(l-2af + d-l] 



q2 qi2 



1 . HqVq") 



(41) 



Writing = 2 — e, expanding around e = 0, replacing 1/e with In/iMS and integrating over a 
we obtain 



d'zi2 K^^^{Xo,Xi]Zi,Z2) 



2Nf f d'^q d^q' 



-iq-{z~xo)+iq' ■(z—xi) 



3tt J (27r)2(27r)^ 



X 



q ■ Q 

q2qi2 



q' In 



q'' In 



q2 _ qi2 



q2 _ qi2 



The details of integrations over q and q' are shown in Appendix ^ The result reads 



d^Zi2K}l^'^{xQ,Xi]Zi,Z2) 



N 



f Z — Xq Z — Xi 



In- 



4e-i-2^ 



2 ' 



67r3|z-a;oP \z - Xi\'^ B?{xq,Xi] z) 
where we have defined a transverse coordinate scale R{xq, xi, z) such that 

In R'^{xo,Xi; z) 



(42) 



(43) 



\z - Xo 


2 In [\z — xi 




■]- 


Z — Xi 


2 In [\z - Xq 






z 


- Xq 


2 _ 


Z — Xi 


2 



or, equivalently, 

R^{Xq,Xi]z) = \z — Xo\\z — Xi\ 



\z 


- iCoP \z - 


- ccip ln(|z - xo\^/ 






- Xo) ■ [z 


-£Ci) \z-Xo\^- 


2 — Xi\'^ 








1 


Xi\ 
Xo\ 




~ (■=-a;o)-(i-!ri) 


2-ii!0l^-|2-a:il^ 



(44) 



(45) 
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Employing Eqs. P^j) and in Eq. we obtain 

Nf Z — Xq Z — Xi 



Q Tt'^ \z — Xq\'^ |, 

4e-f-2^ 



X, 



X < In 



R^{xo,Xi;z) jj:. 



In 



MS 



4e-|-2^ 

2 — Xo\^ 



In 



MS 



4e-t-2^ 

*MS 



(46) 



This is the contribution of the diagrams C in Fig. Q and C in Fig. |21 (along with their mirror 
reflections) to the NLO JIMWLK kernel. To obtain the corresponding contribution to the NLO 
BK kernel one again should use the following formula 



K^'^''ixo,x,;z) 



m,?i=0 



Finally one may substitute the scale R{xq, Xi, z) from Eq. 
obtain 



(47) 

explicitly into Eq. (jlUj) to 



}C^'^''{xo,x,-z] 



N 



f 



6 TT^ 



Xi 



27 



\z - ajol 1^ i^^-^..i2„2 



\z - Xq\ 

|2 



\z - XiY Inr— — 



\Z — Xi\ 



MS 



Xi\ 



ln(|z - Xq 


V 


Z-Xi 






Z -Xo 


2 




Z — Xi 


2 



(48) 



The result in Eq. fl48|) agrees with the NLO correction extracted from the calculation performed 
in [37] where the dispersion method was used in calculating the virtual part of the evolution 
kernel to determine the scale of the running coupling for small- a; evolution. 



3 Ultraviolet Subtraction and Scheme Dependence 
3.1 Subtraction for the JIMWLK Equation 

To understand how the diagrams calculated above translate into corrections to the JIMWLK 
equation, let us recall how the JIMWLK Hamiltonian relates to the leading order diagrams. 
The leading order JIMWLK Hamiltonian is a sum of real and virtual contributions defined 

by 

K'ai[f/] :=Y / d-'xd'ySz /C^«(x,y;;2) ^/^^(zV^zVj, + zV^zV^) ^^g^ 
n\?,,UU] j d'xd'yd'z K^'^ix^y-z) (^V^zV^^ + ^V^zV^) . 

Alternatively we will employ a notation in which an integration convention over repeated trans- 
verse coordinates is implied and write more compactly 

V^^[U] = ^/C^o , Ki^l^l + ^V^zVy + (zV^zV^ + zV^zV^)] . (50) 
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Integration conventions will be implied throughout when we employ subscripts to list the trans- 
verse arguments of the kernels. 

In the above, and are functional derivatives with respect to the path ordered ex- 
ponentials (corresponding to the left and right-invariant vector fields on the SU{Nc) group) 
defined operationally via 

zVlUy := - f/.rCi , tVlUl := rf/t5g) (51a) 

and 

zVlUy ■.=eUJ^^l , iWlUl := -t/irCi . (51b) 

K}"^ was already given in ((221) • notation here is somewhat different from the usual in that 
we absorb a factor l/vr^ into the leading oder kernel. 

The JIMWLK Hamiltonian determines the Y dependence of expectation values of arbitrary 
functionals 0\U] of Wilson lines Ux 

{0[U]){Y):= I D[U]0[U]Zy[U] (52) 

via the Y dependence of the functional weight Zy[U]. The evolution equation for Zy[U] is 
known as the JIMWLK equation: 

dy Zy[U] = -n^^[U] Zy[U] . (53) 

The leading order JIMWLK Hamiltonian in Eq. (j^U)) is constructed such that it adds the 
leading order real and virtual corrections to, say, an interacting qq pair, represented by its 
Wilson line bilinear Uxn ® Ul^ : 




-|- real and virtual self-energy-like terms. (54) 
Taking a trace of Eq. and normalizing by the number of colors turns the above into the 

/ tvU \ 

(unfactorized) right hand side of the BK equation for S{xo,Xi,Y) := ( — 7^"^/ (X) shown in 

trf/ U'^ 

Eq. (|7S)) below. — ^f— ^ is but the most generic of the operators 0[U] referred to in (f^ . 

As in the BK case, real-virtual cancellation and thus UV finiteness follow from the appear- 
ance of the same kernel /C'"'^ in both real and virtual contributions. This ensures that the limits 
z X and z ^ y cancel between the two terms under the integral. Probability conservation 
at leading order manifests itself more globally in the absence of interaction with the target, i.e.. 
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in the limit f/ — ^ 1: There ensures that 'HviS;uai ~^ ~'^rcai that there is no evolution 
without interaction with the target. 

7^^° may in fact be used to act on any tensor product of quarks, antiquarks and gluons 
in a projectile's wavefunction that interact with the target via corresponding Wilson lines to 
produce a sum of leading order as ln(l/xBj) corrections to this eikonal interaction. This is the 
technical mechanism by which the JIMWLK equation translates into the Balitsky hierarchy. 

Note the efficiency with which the JIMWLK Hamiltonian encodes the contributions: due 
to the symmetry properties of the kernels, the self-energy-hke diagrams arise from the same 
terms that create the exchange diagrams. The only distinctions left in the Hamiltonian are: 

• The order of the vertices w.r.t. the target interaction, i.e. the Wilson lines. This is 
encoded in the use of the V and V. 

• The interaction (or lack thereof) encoded in presence or absence of an adjoint Wilson line 
at the transverse position at which the newly created gluon interacts with the target, as 
shown in the second and third lines of Eq. ()49|). 

This pattern extends itself to the NLO contributions studied here. Only the variants of dia- 
gram A differ slightly in structure from the contributions already encountered at leading order: 
they depend on two new transverse coordinates and contain a factor 2tT{t^Uzit"'Ul^) instead of 
the U^^ of the real emissions at leading order. This leads to the following correspondence of 
diagrams and terms in the NLO corrections in the Hamiltonian 

ln(l/xBj) U^o ® = ln(l/xBj) J (fx(fyd''zi(fz2 /CfLO(a;, y; Zi, z^) 




(55a) 

All other corrections have the same U structure already encountered at the leading order. We 
have corrections to real emission 

ln(l/xBj)7^f° U^,®Ul=a^ ln(l/xBj) J Sxd^y(fzlC^^''{x,y-z) 

X Uf (zV^zVj, + iWl^D ® Ul = + + ... (55b) 
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and virtual terms 



ln(l/xBj)H^^° U^o^Ul = -a^ ln(l/xBj) J (fx(fy(fz IC^'^''{x,y; z) 

i^l^l + ^ v:^ V^) f/.„ ® ul 



X 




+ . . . 



(55c) 



The minus sign in the last term is due to the different V structures in real and virtual terms 
and is important for the real-virtual cancellations. The dots represent both symmetrization in 
external coordinates Xq and Xi as well as the inclusion of "self energy like diagrams" in which 
the gluon line connects back to the quark (or antiquark) it originates from. 

We group the contributions accordingly (again employing an integration convention for all 
repeated transverse coordinates x, y, z, zi) 



new: 
real: 
virtual: 



NLO 



HI 



-i/LO I -T/NLO 
'T-real T~ ''-2 



a 



M r-NLO 



2 '^l x,y;zi,Z2 



virtual 



LO 
LO 



2tr(t''f/.,rf/tj(zV:zVj, + zV:zV; 



(56) 
(57) 
(58) 



and observe also at this order, that, due to probability conservation as expressed by Eq. ()H6|1 . the 
limit f/ — > 1 leads to a cancellation of the sum of all these contributions. We note: probability 
conservation connects all of the above contributions. 

The above separation of terms is quite unsatisfactory also if we wish to extract the running 
coupling contributions to the leading order Hamiltonian. Two complementary issues emerge: 

1. Any running coupling correction should come as a uniform modification in both real and 
virtual terms of the leading order kernel, replacement 



a^/CLO(a;, y; z) ^ a^/C^"(^, y- z) + «^/C^^"(x, y- z 



-LO 



2 irNLO/ 



(59) 



with a yet unspecified NLO kernel. This is required if an all orders resummation of quark 
bubbles is to take the form 



a^/CLO(a;, y- z) ^ a,(/(a;, y; z))K}^''{x, y; z) (60) 

inside the coordinate integrals of the JIMWLK Hamiltonian and if the pattern of real 
virtual cancellation (and thus probability conservation) be maintained beyond the leading 
order. The sum of real and virtual contributions in the above is not of this form; there is 
no common NLO kernel in both terms. Not even the divergent contributions (traceable 
by the //-dependence of the transverse logarithms) in Eqs. (jHT|) and (jlHj) coincide. This 
is related to the second issue: 
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2. The new term, Eq. ()55ap . contains UV divergent contributions where I212I, the separation 
of quark and antiquark, reaches the UV cutoff. To extract the UV divergence, which is 
driven by scales much larger than the saturation scale Qg, the ^/-dependent part of the 
quark loop, the factor 2tT{t'^Uzit°-Ul^), may be expanded in Z12 around some fixed base 
point z 



2ti{t'u,,eul) = u\ 

+ [{z - zi)dz 



Z2)9,j2tr(t^[/,,rf/t ) 



+ 



Zl.2=Z 



(61) 



so that to leading order in this expansion (i.e., keeping the term only) the resulting 
[/-dependence of ()55a|) takes a form similar to that in ()57|) . The integral over Z12 may then 
be carried out and its divergence exposed as illustrated in Fig. 01 [For this to be sufficient 




UV-divergcnt 



Figure 3: Separating UV-finite and UV-divergent parts of Fig. QA 

it is of course mandatory that only the leading order in this Taylor expansion contains 
a UV divergence.] This divergent contribution must patch up the mismatch between the 
real and virtual terms discussed previously. While the divergence is independent of the 
choice of base point, the finite terms associated with the separation shown in Fig. El will 
depend on this choice. This will lead to a scheme dependence to be discussed below. 

For the JIMWLK Hamiltonian, we are thus led to consider a term of the form 



a. 



Jn(l/xBj) jdh^dh,2 }C'^''''ix,y-z,,Z2)UfitVltV'y + tVlt\/'y) U^^^U^ 




(62) 



which carries the UV divergence of (j55ap . 

By subtracting this contribution from ()56p and adding it to (|57p. we shift the UV divergence 
from a genuinely and physically new contribution in which a distinguishable, well separated qq 
pair interacts with the target, to the contribution that is not distinguishable from the single 
interacting gluon already present at leading order. While the logarithmically UV divergent term 
is uniquely defined, the finite scale dependent terms under the logarithm are not constrained. 
This is the origin of our scheme dependence. 
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To be explicit, we make use of ()38p for our choice of z and define what we will call the 
subtraction term in the gluon scheme* 

where the kernel in the subtraction is calculated with z placed at the gluon position z: 

tf^'^ix.y-z) = I d^z.^lC'^^^^ix^y-z.^z,). (64) 

The exphcit form of the right hand side of Eq. (jMj) was aheady obtained in our calculation of 
the fully virtual corrections in Sect. 12.31 with the answer given by Eq. (j43|l . We use it first to 
define a genuinely UV finite qq contribution of the form 

-^l^i'£y;.,,., 2tr(t"f/.,t^f/tj(zV:zV^ + zV^zV^) - ^/Cf^Ji^, U:\^l^Vl + zV^zVj,) . 

(65) 

This contribution only is of interest if we wish to go beyond the inclusion of running coupling 
corrections to include genuine NLO contributions. While such a calculation would be interesting 
and important, it remains beyond the scope of this paper. Here we note that the term in Eq. ()65p 
is UV- finite and vanishes by itself in the no- interaction limit oi U —>■ 1: This term no longer 
mixes with the remaining contributions under probability conservation and -contrary to the 
unsubtracted contributions- we may neglect if we are only interested in the the scale of the 
running coupling. 

The remaining contributions now assemble directly into a form that fulfills all requirements 
of a running coupling contribution. Using (jSHI), we find that adding (jMjl to (jHZj) leaves us with 
identical kernels both for real and virtual contributions 

a, /CLO(x, y; z) - /Cf ^(a;, y; z) = a, l&^'ix, y; z) + al (^f^°(a;, y; z) + /C^O(a^, y, z)^ . 

(66) 

The leading Nj contributions to the running coupling corrections for the JIMWLK Hamiltonian 
take a from analogous to (j^Uj) 

Uf{^Vl^Vl + zV^zV^) + (zV^zV^ + zV^zV^)" . 

(67) 

The equations in the Balitsky hierarchy created with this operator are finite and unitary for 
fixed projectile configurations. The subtraction fully decouples conformal contributions ()65p 
and non-conformal contributions ()67|) up to the order making it feasible to discuss running 
coupling corrections independently of the conformal contributions. 

*The concept of an explicit UV subtraction was first introduced by Balitsky in the calculation of transverse 
coordinate space version of NLO BFKL in [49]. We thank Ian Balitsky for communicating it to us in private. 



a. 



LO 



a. 



LO 



I i-LO 
"T '^2 x,y;z 
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3.2 Subtraction for the BK Equation 

The UV subtraction described above for JIMWLK evolution equation can be translated to the 
BK framework [17-19] simply by repeating the steps that allow to identify BK as a limiting 
case of JIMWLK at the leading order. Here we will instead formulate the argument again 
entirely within the BK framework to provide a self contained discussion. We begin by writing 
the standard LO BK evolution equation for the forward amplitude of a quark dipole scattering 
on a nucleus 

Nix,, x„Y)^l-^^ (tr [U^, f/t J ) (F), (68) 

where f/'s are from Eq. |TT| . the transverse coordinates of the quark and the anti-quark are Xq 
and Xi, and the dipole's rapidity is Y. The LO BK equation reads 

dN{xoXi,Y) ^ I ^^^^^^ ^ ^^^^^ _ ^^^^^ 

-N{xo,X2,Y)N{x2,x^,Y)], (69) 

where Xmn = \Xm — Xn\ and the large- A''^ limit is assumed. Using the LO JIMWLK kernel from 
Eq. (|32p we can define the LO dipole kernel by 

1 

K'^''{xo,x,;z) =Cf (-l)'"+"/CL°(a^™,^n;;z). (70) 

m,n=0 

The dipole kernel (j7U|) sums up the same diagrams as shown in Eq. (j54j) for the LO JIMWLK 
Hamiltonian [22]. Using Eq. ()70p we rewrite Eq. (jU^ as 



dN{x,,x,,Y) 
dY = 



j d^X2K^°ixo,xr,X2) [N{xo,X2,Y) + N{x2,Xi,Y) - N{xo,xi,Y) 



-N{xo,X2,Y)N{x2,x^,Y)]. 

(71) 

For the purpose of performing the UV subtraction, it is more convenient to rewrite Eq. (|7T|l in 
terms of the S-matrix 

S{xo, xi, Y) = ^ (tr [f/,„ f/t J ) (F) = 1 - N{xo, x^, Y) (72) 

obtaining 
dS{xo,Xi,Y) 



dY ~ "'^ 



j (fx2 i^^°(a;o, xi- X2) {Six,, X2, Y) 5(0=2, x^,Y) - Six,, x,, Y)] . (73) 



Now we are ready to include the NLO corrections calculated in Section |21 Adding the 
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diagrams in Figs. ^ and |H to the LO dipole kernel yields the following evolution equation 



dS{xo,xi,Y) 
dY 




{xq, xi, X2) [S{xo, X2, Y) S{x2, xi, Y) - S{xo, xi, Y)] 




{xq, Xi, Zi, Z2) S{xq, zi, Y) S{z2, Xi, Y) 



I 



{xo, Xi; X2) S{xo, X2, Y) S{X2, Xi, Y) 



{xq,Xi;x2) S{xq,Xi,Y). (74) 



Similar to the JIMWLK case we notice that while kernels and Kf^'^ appear as higher 

order corrections to the leading order BK kernel K^'^ having the same transverse coordinate 
dependence, the kernel K^^'^ stands out. It includes integrals over two transverse vectors, Zi 
and Z2i instead of one. Since the shape of the leading kernel K^^ in not preserved in K^^^ 
and as we are looking for running coupling corrections to K^'^ , one may naively discard K^^^ 
as not giving any running coupling contribution. However, before we embark on extracting the 
running coupling corrections, let us formulate general rules for such corrections. Similar to the 
JIMWLK case we require the following: 

• Unitarity: as Eq. (fTSj) gives an explicitly unitary solution for S, i.e., as rapidity Y 

00 then S* — s> 0, we require the running coupling corrections to preserve this unitarity 
property. This requirement is satisfied as long as the right hand side of the evolution 
equation has only terms containing powers of S. 

• No interaction — no evolution condition: we require that in the absence of interaction 
the right hand side of the resulting evolution equation should become when 5" = 1 is 
inserted there. This condition is easily satisfied in the standard Feynman perturbation 
theory where the non-interacting graphs are zero. In the light cone perturbation theory 
(LCPT) the non-interacting diagrams are not zero, which allows us to define and calculate 
light cone wave functions. Because of that it is a little harder to show in LCPT that in the 
absence of interactions all diagrams for the amplitude (not to be confused with the wave 
function) cancel. For instance, the no interaction — no evolution condition is satisfied by 
Eq. ()74|) : if we put 5* = 1 on its right hand side we will get zero due to the condition in 
Eq. (j2Hl)- We want this property to be preserved after running coupling corrections are 
included. 

From the above conditions one can see that simply discarding Kf^'~' from the right hand side 
of Eq. (j74j) would not work: while the equation obtained this way would satisfy the unitarity 
condition, it would not satisfy the second condition stated above, since, for 5* = 1 we will not 
get zero on the right hand side anymore. What strengthens the case for keeping a part of Kf^*^ 
is that it contains a UV divergence, as can be seen from Eq. (pSj) . which may contribute to 
the running of the coupling constant. Kernels ii'^'"^ and Kf^'^ also contain UV divergences, 
which need to be canceled by the divergence in Kf^*^ as follows from Eq. (j38|) . Without Kf^*^ 
the right hand side of Eq. (|74p would become infinite. Therefore, to keep the right hand side 
of the resulting evolution equation finite, and in order to satisfy the second one of the above 
conditions, we propose the subtraction illustrated in Fig. EJ 
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UV-finite 



UV-divergent 



Figure 4: Separating UV-finite and UV-divergent parts of the kernel in the NLO BK 

evolution. The ovals denote color dipoles. 

Formally we write 

j d^zi d^Z2 i^r"(a;o, xr, Zi, 22) S{x^, Zi, Y) S{z2, x^, Y) 

= j dhd'zuK^''''ixo,x^,z^,Z2) [Sixo,z^,Y) Siz2,x,,Y) - Sixo,z,Y) Siz,x^,Y) 

+ Sixo,zX)Siz,x^X)], (75) 

where z is the position of the virtual gluon in Fig. defined in Eq. Q.^ Now the first two 
terms in the square brackets on the right hand side of Eq. ()75p give a UV-finite result, as 
shown in Fig. HI which goes to zero both for S* = and S = 1. These terms combined do not 
have a UV- divergence and do not contribute to the running coupling constant. They give a 
non-running coupling NLO BK evolution piece and we will discard them here. The last term 
on the right hand side of Eq. ()75|) we will keep. Similar to the JIMWLK case we define the 
subtraction kernel by 

Kf'^''{xo,x^,z) = J d'z,,Kf'^''{xo,x^,z,,Z2). (76) 

An explicit form of Kf^'-'{xQ, Xi; z) can be found using Eqs. ()43p and ()16p . With the help of 
the definition in Eq. ()76|) we rewrite the last term in Eq. (fTSj) as 

j d'x2 i^r°(^o, ^1; X2) S{xo, X2, Y) Six2, Y). (77) 

t Indeed as z from Eq. depends on the longitudinal momentum fraction of the quark a, switching from 
Zi and Z2 to Zyi and z implies a change in the a-integral in Kf^^ . However, since above we have used 212 and 
z vectors everywhere, no changes apply to our earlier results. 
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Keeping only this term in kernel Kf^'^ modifies Eq. (f7^ to give 



dS{xo,Xi,Y) 
dY 



I 




(a?o, xi, X2) [S{xo, X2, Y) S{x2, xi, Y) - S{xo, xi, Y)] 





I 
I 



(f'X2 -ft'f^°(a;o, xi] X2) S{xo, X2, Y) S{x2, xi, Y) 



{xq, Xi, X2) S{xo, X2, Y) S{x2, xi, Y) 




I 



{xq,Xi;X2) S{xq,XiX)- (78) 



Using Eq. (jHHj) we rewrite Eq. (fTHj) as 



dY 



I 



d^X2 K^'^ixo, xi; X2) + aj i^f ^°(a;o, Xi, X2) + al K2^'^{xo, Xi, X2) 

X [S{xo, X2, Y) S{x2, xi, Y) - S{xo, Xi, Y)] . 

(79) 



Eq. (ffn|) obeys both of the conditions stated above: its right hand side is zero at both 5 = and 
5 = 1. Moreover the right hand side of Eq. (ffn|l is UV finite, which is essential for obtaining a 
meaningful result. The kernels Kf^'^ and K2^'~' both look like corrections to the LO kernel. In 
the following, when we study fermion bubble insertions to all orders, we will use the format of 
Eq. (fTHj) to systematically include their contributions into the running of the coupling constant. 

The choice of subtracting and adding S{xq, z, Y) S{z, Xi, Y) depending on gluon's position 
z in Eq. (j75p is indeed quite arbitrary. For instance, one can use Zi or 22 (or any other linear 
combination of the two vectors Zi and Z2) in place of z} We can not find any argument or 
criterion which would prefer one choice of the "subtraction point" over the other. We choose 
z as our "subtraction point" since it appears to be convenient and goes along the lines of cal- 
culating A'!^^*-* in Eq. This choice appears to also be preferred by the dispersive method 
of calculating the running coupling correction to small-x evolution used in [37]. Indeed the 
uncertainty in selecting the "subtraction point" does not affect our ability to extract the UV 
divergent part of Kf^'^. However, it may change the scale R under the logarithm in Eq. ()43p . 
resulting in a different scale for the running coupling constant. We believe that the modifica- 
tion of the running coupling scale due to varying the "subtraction point" will be numerically 
insignificant: however, a detailed study of this question is left for further investigations. Here 
we will refer to this dependence of the running coupling scale on the "subtraction point" as of 
some sort of a scheme dependence for the running coupling constant. 



^We thank Ian Balitsky for helping us to reach this conclusion. 
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4 Resummation of Bubbles to All Orders: Setting the 
Scale for the Running Coupling Constant 



Now we are ready to resum all powers of as Nj corrections in the JIMWLK and BK evolution 
kernels. To accomplish that one has to insert infinite chains of gluon bubbles onto the gluon 
lines in Figs. ^ and |21 An example of corresponding higher-order diagrams is shown in Figs. El 
and El 




ABC 



Figure 5: Diagrams giving the higher order a^Nf corrections to the kernels of JIMWLK and 
BK small-x evolution equations. To get the all-order UgNf contribution one has to sum an 
infinite series of quark bubble insertions. 




A' C 

Figure 6: Diagrams giving the higher order OgNf corrections to the kernels of JIMWLK and BK 
small-x evolution equations containing instantaneous gluon lines. Again, to get the all-order 
agNf contribution one has to sum an infinite series of quark bubble insertions. 

An explicit calculation using the rules of light-cone perturbation theory [41,42] shows that 
inserting all-order quark bubbles on the gluon lines generates geometric series in momentum 
space. Before calculating the diagrams in Figs. El and El we remember that, as was discussed 
above, in order to find the running coupling correction, instead of the diagram in Figs. EK and 
EK' we should consider the "subtraction" diagrams A and A' shown in Fig. [7| 
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A A' 

Figure 7: The "subtraction" diagrams which should be considered in place of the diagrams in 
Figs. and\^A' for the running coupling scale calculations. 



By an explicit calculation, similar to the calculation of a single bubble insertion which led 
to Eq. ()28|). one can show that the contribution of the "dressed" subtraction diagrams A and 
A' in Fig. [7|to the JIMWLK kernel reads 



X 




X -. — w T^' (80) 



where we have also replaced all factors of Nf by —Qir (32- Similarly, the contribution of the 
"dressed" diagram in Fig. El B is obtained by iterating the quark bubbles from Eq. ()28p on both 
sides of the cut. The series of bubbles on each side of the cut generates a geometric series. The 
zeroth-order term in this series is the leading order JIMWLK kernel. The first-order term in 
the series is given by in Eq. (jSH). It is therefore more convenient to write down the sum 

of the LO JIMWLK kernel and the contribution of the diagram in Fig. El B. The result is 



X -. ^. (81) 

MvIS" / \ MvIS" 



The contribution to the virtual part of the JIMWLK or BK kernels is given by the diagrams in 
Figs. andinp'. It can be easily extracted from Eqs. (|Hn|) and (|HT|) using the condition 
which holds to all orders in quark bubbles. The resulting virtual kernel would be equal to the 
real kernel as was explained above and shown in Eqs. (|67p and (|79p. The dispersive method 
of [37], which was used there to calculate Figs. andlHP', can also be used to obtain Eq. (jHTj) . 
and, by employing the probability conservation condition (j^Hj) . to recover Eq. (jHUj) as well. 
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Before proceeding to evaluate the kernels in Eqs. ()8U|) and ()81|) let us first analyze their 
dependence on the UV cutoff /i^s- This is instructive because the cutoff is indeed a constant 
and therefore /i^^-dependence does not get modified by the Fourier transform. Hence the ftj^- 
dependence of the kernel is the same in transverse coordinate and momentum spaces. Keeping 
only a^, P2 and /i^jg we thus write 

«^/C® oc , (82) 

and 

From Eq. (j83p we immediately see that the sum of the LO kernel and the diagram in Fig. 03 does 
not give us a renormalizable quantity, as it can not be expressed in terms of the renormalized 
coupling constant. It lacks a power of to give us a square of the physical coupling a^. Now 
the need of extracting the UV divergence from the graph in Fig. becomes manifest. Adding 
the "subtraction" term (|H11) to (|HS1) yields 



r-LO , 2 r- , 2 (1 - a^/32ln/iy a, a, 

a,/C +a,/C@ + .,/C(r)Oc _ ^^^^^.^ oc --, (84) 

which is indeed renormalizable as it can be expressed in terms of the physical coupling a^. (We 
have left a factor of both in the numerator and in the denominator of Eq. (jS^ on purpose 
to underline the fact that the arguments of all three couplings, which we did not keep, may be 
different, which would prohibit the cancellation.) 

The sum of the kernels in Eqs. ()8()j) and (jHT|l gives the all-order in UsNj (or as/52) contri- 
bution to the running-coupling part of the real JIMWLK kernel: 

a^,}C,c{xQ,Xi]z) = a^}C^'^{xo,Xi;z) + al}CQixo,Xi;z) + al}C(r){xo,Xi;z) = 

-5/3" 



rl^n rl'^n' n n' ft^ 1 + a^/32 In ■ 



-p. 

MS 



(2^)' (2^)' (l + a,i5. In ^) (l + a,/?2 In ^ 

\ ^MS / V MvIS" 



(85) 



where we have defined a momentum scale Q by 



(Indeed is independent of be easily seen from Eq. (jHUj) .) As the one-loop running 

coupling constant is defined in the MS scheme by 

rv 

o^sm = ' (87) 

1 + a^/32 In ^ 
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Eq. ()85|) can be rewritten as 

y (27r)2(27r)2 ^2^/2 a,(g2e-5/3) " 

Eq. (j88p is the first of tfie two main results of our paper. It gives the JIMWLK kernel with 
the running couphng corrections included in transverse momentum space. Remarkably, the 
corrections come in as a "triumvirate" of the couplings,^ instead of a single coupling constant 
with some momentum scale! Despite the surprising form this result is in full agreement with the 
expressions found in [37]. To facilitate comparison, we give a detailed translation in appendix. O 
This provides for an interesting mechanism to reduce the above result to a simpler underlying 
structure expected for the purely virtual contributions of diagrams Fig. El C andlHlC': There 
the z integral may be performed and, in the absence of interaction with the target, this sets 
q' = q. Then the sum of the integrands of transverse and longitudinal contributions reduces to 



ia.(g2e-t) (89) 



which clearly corresponds to the exchange of a dressed noninteracting gluon. This is the 
counterpart of a cancellation that appears in [37] under the same premise. Note that due to 
the virtual nature of the diagrams in Figs. El C and El C the running coupling corrections can 
enter only with the sole available transverse momentum scale q^ in its argument, as shown in 
Eq. ()89|) . Therefore, while different subtraction procedures may yield different expressions for 
the scale Q, as compared to Eq. (jHEl), all of these alternatives must lead to expressions for 
that approach q^ in the limit when q' = q. This would reduce the "triumvirate" of couplings 
from Eq. (jHH|) to the single coupling shown in Eq. ()89p. as expected for virtual diagrams here. 

The most intriguing property, however, is that the "triumvirate" structure of ()88|) solves 
the puzzle of how to successfully perform a BLM scale setting: In [37] it was observed that 
an attempt to perform a BLM scale setting in a perturbative formulation with a single Borel 
parameter (corresponding to an approximation in terms of a single geometric series in our 
present language) would not lead to a successful resummation of the dominant contribution. 
From the present perspective this would correspond to the attempt to use the joint leading 
order expansion of all three couplings in (jHH|) to determine a single scale (instead of the three 
separate ones of ()88|) ) that would give a good approximation to the triumvirate in terms of 
a single geometric series. Any such attempt would necessarily entail an artificial all orders 
iteration of the "numerator logarithms" encoded in the "denominator coupling" as{Q'^e~3) 
that is clearly absent in the underlying expression. This is the source of the spurious and 
divergent higher inverse powers of q- q' or {xq — z) ■ [xi — z) encountered in a naive attempt of 
deriving a BLM approximation in [37]. The "triumvirate" structure will allow for a successful 
and transparent BLM approximation of the full coordinate result given in [37]. 

''We note that a similar structure containing three couphng constants has been obtained independently by 
Balitsky [49] . Our difference from [49] appears to be due to a different choice of the "subtraction point" . 
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Eq. ()88|) allows one to find the corresponding BK evolution kernel with the running coupling 
corrections included by using 

1 

K,,{xo,Xi;z) = Cf Yl {-ir^''ICUxm.,Xr,;z). (90) 

m,n=0 



To determine the scales of the three physical couplings in Eq. (jHHjl in transverse coordinate 
space we have to evaluate the integrations in the kernels from Eqs. fl8()|l and (jHTjl . In transverse 
momentum space each chain of bubbles generates a geometric series, as shown in Eq. ()85|1 . 
However, a Fourier-transform of these series is dangerous for several reasons. First and foremost 
the integrals over q and q' in Eqs. fj80|) and (jHT|) include contributions from the Landau poles 
leading to power corrections which are not under perturbative control. The uncertainties due 
to power corrections are estimated in [37] using renormalon techniques [50,51]. Our strategy 
here is to ignore these contributions concentrating on setting the scale of the running coupling 
in transverse coordinate space. Even then our goal is difficult, since, even though leading 
powers of In terms generate a geometric series in transverse coordinate space just like in 
the momentum space, it is not clear whether the transverse coordinate scale under the logarithm 
stays the same in all the terms in the series. In that sense, setting the running coupling scale 
in transverse coordinate space will only be an approximation of the more exact Eq. (jHH|) in the 
sense of a BLM scale setting. 

To study running coupling corrections in transverse coordinate space we begin by evaluating 
the Fourier transforms in Eq. ()81|) . Similar to Appendix El we perform the angular integrals 
first to write 



a^}&'^{xo,Xi;z) + allC@{xo,Xi;z^- ^ 



an z — Xq z — X] 



X / dqdq' Ji{q\z — Xq\) Ji{q' \z — Xi\) , . , 

/ n 1 <72p-5/3\ „ , „/2„-5/3 

J l + a^/52ln^-^ l + a^/^sln^-^ 



(91) 



MS / V ^MS 

with q = \q\ and q' = \q'\. Since our goal is to find the scale of the strong coupling constant 
ignoring the power corrections we can expand the denominators of Eq. ()91|) into geometric 
series obtaining 



a,]C'^''{xo,x,;z) + al}C@ixo,x^;z) = ^ . J- {-a.^.T^^ 

n \z — Xo\ \z — Xi\ 



X dqdq' Ji{q\z-Xo\)Ji{q'\z-Xi\)\Yi''^^ In™ ^-^ . (92) 

Rewriting the powers of the logarithms in Eq. (jH21) in terms of derivatives yields 

oo 

a^l^'^^x.^xr-z) + alK:@{x,,xv,z) = ^ . J- ^-a.P^T^^ 

\Z — Xo\ \Z — Xi\ 

d"" d"^ / g^e-^/M / g'2 g-5/3 



d"' d^ f 

^dF'dX^ j '^'i'^l' JMz-XG\)Ji{q'\z-Xi\ 







(93) 
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Performing the q- and g'-integrals gives 

oo 



n+m 



n,m.=0 



4e-^/3 \ r(l + A) 



^-^ol^Ms/ r(i-A) 



A=0 



4e-5/3 



r(i + y) 



A'=0 

(94) 



Differentiating with respect to A and A' we write out the first few terms in the resulting series 
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(95) 



One can see that the geometric series structure appears to hold up to the cubic terms in either 
one of the logarithms. In the sense of a BLM-type approach [43] we approximate the expressions 
in each of the curly brackets by a geometric series, obtaining 



a^}&'^{xo,Xi; z)+alK:@{xo,Xi; z) ^ }&'^{xo,Xi; z) 
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(96) 



Evaluation of the Fourier transforms in Eq. (j8Up is performed along similar lines in Appendix 
IbI The result reads 



/ 4e"l"2^ 

^1^®{xq,xi]z) ^ K}'^{xQ,Xi]z)a^(32 In — 
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(97) 



with the scale R given by Eq. pHjl . 

Finally, adding Eqs. and (|U7j) yields 



a^/Crc(xo,a;i; z) ^ /C^°(xo, a;i; z) 



4„-5/3-27 \ / . 5/3-27 



4g-5/3-27 
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with the scale R{xq,Xi; z) given by Eq. (jl^j) and the couphngs calculated in the MS scheme. 
Eq. (|98p is the second of the two main results of our paper. It gives the JIMWLK kernel with 
the running coupling constant. It is very interesting that the running coupling corrections come 
in not through a scale of a single running coupling as one would naively expect, but in the 
form of a "triumvirate" of the running couplings shown in Eq. (j^Hj) ! Eq. ()98|) can be used to 
construct BK kernel with the running coupling constant by employing Eq. (jMIj) . 



5 Conclusions 

To conclude let us state the main results of this work once again. By tracking the powers 
of as Nf we have included running coupling corrections into the JIMWLK and BK evolution 
equations. Our all orders result agrees with the expressions derived with the dispersive method 
in [37], but renders the result in terms of a "triumvirate" of couplings already in the momen- 
tum space expressions (jHH|l . This allows us to give a concise, accurate BLM approximation of 
the perturbative sum in coordinate space in terms of a corresponding coordinate space "tri- 
umvirate" shown in Eq. ()98|) . Our procedure includes one uncertainty, related to choosing the 
"subtraction point", as was discussed in Sect. El This ambiguity is akin to scheme- dependence 
of the running coupling constant and we believe that the final result does not depend on the 
choice of the "subtraction point" in a very crucial way. We picked the subtraction point to be 
at the transverse coordinate of the virtual gluon in Fig. ^A.. 

Our result for the JIMWLK Hamiltonian with the running coupling constant is 

(99) 

with integrations over x, y and z implied. All the running couplings should be calculated in 
the MS scheme. 

To obtain the BK evolution equation with the running coupling constant one needs to sum 
the kernel in Eq. (j98p over all possible connections of the gluon to the quark and the anti-quark 
lines, as formally shown in Eq. (jUU)) . (Alternatively to derive the BK equation one can apply the 
JIMWLK Hamiltonian from Eq. ()99|) to a correlator of two Wilson lines and take the large- iVc 
limit.) A straightforward calculation shows that 

lim R'^{xo, Xi; z) = \z — Xq]'^ and lim R'^{xo,Xi]z) = \z — X\^ , (100) 

such that the BK evolution equation with the running coupling corrections is 
9iV(a;o,xi,F) Cp 



dY vr 
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, 4e-^/^-^^ ^ 1 ^ «^ l^S— j l^j— ; ■ ^21 , /4e-5/3-2n i 
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20 / -^20 — ■ r •^20-^21 \ -^21 / -^21 



X [N{xo, X2, Y) + N{X2, xi, Y) - N{xo, x,,Y) - N{xo, X2, Y) N{x2, x,, Y)] , (101) 
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where is given by Eq. ()45p. 

It is interesting to explore the limits of Eq. (jlUlj) . We will refer to the original dipole 01 as 
the "parent" dipole, while the dipoles 20 and 21 generated in one step of the evolution will be 
called "daughter" dipoles. First of all, when both produced dipoles are comparable and much 
larger than the "parent" dipole, ~ X21 ^ xoi, the argument of the coupling constant in all 
three terms in Eq. pOlj) would be given by the "daughter" dipole sizes 0:20 ~ 2:21 • However, 
such large dipole sizes should be cut off by the inverse saturation scale 1/Qs, which implies 
that the scale for the coupling constant in the IR region of phase space would be given by Qs 
keeping the coupling small and the physics perturbative. In the other interesting limit when 
one of the "daughter" dipoles is much smaller than the other one, X20 ^ 3^21 ~ a^oi; a simple 
calculation shows that R^{xq^ Xi] X2) ~ and the BK kernel in Eq. pOlj) becomes 

\ "''20 / "''20 V "''21 / "''20 "''21 V "''21 / "''21 

Two out of three terms in the kernel have the scale of the coupling given by the larger dipole 
size, naively making the evolution "less perturbative". However, in the X20 ^ a^2i limit it is 
the first term which dominates Eq. ()102|1 : that term has the running coupling scale given by 
the size of the smaller dipole, making the physics perturbative! 
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A Evaluating the Fourier transforms in Eq. (1421) 



Here we first calculate the following integral coming from the first (transverse) term in Eq. 



^''^'''^^ - J j2nYj2^^' ' ^^^^ 

where z and z' are some transverse coordinate vectors. Replacing q ^ idz and q' — > —i dz' in 
the numerator of the first ratio in the integrand we can integrate over the angles of q and q' 
obtaining 

1 o o M?^?' w,... ^j-^ (,1^S= 3 



I^(z,z') = -—dz-dz> ^^J,[qz)Uq'z')—^^ L ^ ^ L (A2) 

[2ti)^ J q q q'^ — q^ 
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with q = \q\, q' = \q'\, z = \z\ and z' = \z'\. Bringing the transverse gradients back into the 
integrand yields 



It{z,z') = I dqdq'Mqz)Mq'z') (A3) 

{2n)'^ zz' J q^ — q'^ 



To perform q and q' integrations we first rewrite the fraction in the integrand of ()A3|) as 

^2 72 ^ = 1^^-^ + ?"^^ (A4) 

q2 - q,'2 ^ q2 _ qi2 

and then use the following integral identity 
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/ g2(i_^)+g/2^ q2 _ q,2 


to replace the last term in Eq. ()A4|) .^ Eq. ()A3|) becomes 

cx> ( 1 

It{z,z') = / dqdq'Mqz)Mq'z') hn + g'^ / dp 

[ZTT) ZZ J I fXj^ J 
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q^{l-P) + q'^P 

(A6) 



Now the q and q' integrations can be performed using standard formulas for the integrals of 
Bessel functions yielding 

1 z ■ z' \ 4 e't"^"^ /" 1 
It{z,z') = ——^^ Ihi di3 y (A7) 

Finally the /9-integral in Eq. ()A7jl can be done using Eq. ()A5|1 giving 



1 _ _/ -5 m ,2 2 2 m 2 ,,2 — 

- ^ 

Now let us use the same technique to perform Fourier transforms in the second (longitudinal) 
term in Eq. (j42j) . We want to evaluate 

= / e---»-' !^^, (A9) 

J (27r)2 (27r)2 - q'^ 

First we integrate over the angles of q and q' 

00 

h{z,z') = dqd<jqq'Mqz)Mq'z')^^^^. (AlO) 



^We thank Ian Balitsky for pointing out to us the uscfuhicss of this substitution. 
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Now we can use Eq. ()A5|) to integrate over q and q' obtaining 



1 







which, using Eq. ()A5|1 again we can write as 



B Evaluating the Fourier transforms in Eq. ( 1801) 



Here we will try to perform the Fourier transforms in Eq. ()80|). We begin by analyzing the 
transverse part of the kernel, given by the first term in the curly brackets in Eq. (jHUj) . We start 
by performing the angular integrations over the angles of q and q', which, similar to the way 
we arrived at Eq. ()A3|) . yield 

c? K7^{xq,Xi\z) = \ I dqdq' Ji{q\z - Xq\) Ji{(^ \z - Xi\) 

TT^ \Z — Xq\ \Z — Xi\ J 
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with /Cq denoting the transverse part of the kernel. Since our intent is to extract the scale of 
the running coupling in transverse coordinate space ignoring power corrections, we expand the 
denominators in Eq. ()B1|) into geometric series and repeat the steps which led from Eq. f)A3|) 
to Eq. ()A6|) writing 



00 

2 a 00 „ 
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(B2) 

It is hard to perform q and q' integrations for a general term in the series characterized by some 
values of n and m. However, we can try estimating the first correction, i.e., the n = l,m = 
term: 

C(fif32 Z — Xq Z — Xi 
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(The n = 0, m = 1 term will be constructed by replacing a:;o ^ a^i in the result of evaluating 
fjB3|) .) The g'-integral can be easily done in Eq. ()B3|) . along with the g-integral in the first term 
in the brackets, yielding 
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To perform the g-integral we first re-write the remaining logarithm as a derivative 
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Performing the g-integration yields 
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Using the definition of hypergeometric functions we write 
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With the help of Eq. ()B7|) the differentiation with respect to A can be easily carried out in 
Eq. ()B6|) . After integrating over j3 we obtain 
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With the help of "transverse" part of Eq. ()43|). in which we replace Nj —671^2, we derive 
the following expansion: 
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Here the ellipsis include not only the higher order terms in but also the term quadratic in 

with the Xq x^ replacement. It can be shown that the term in the last line of Eq. ()B9|) 
is numerically small compared to the other terms in the series. Dropping that term yields 
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It appears likely that the higher order corrections would continue the geometric series in 
Eq. (|BT0|1 with the same constants under the logarithms. Resummation of such series yields 
the "transverse" part of Eq. (^7}) . 

Now we have to evaluate the longitudinal (instantaneous) part of Eq. (jSOl), given by the last 
term in the curly brackets in that equation: 
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Again the exact integration does not appear possible. Instead we will expand the running 
coupling denominators to the linear order in the logarithms and evaluate the following term 
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which we rewrite using Eq. flA5j) as 
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Performing the g'-integral first yields 
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where we have again replaced the logarithm by a derivative of a power. Integrating over q we 
obtain 
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Using the expansion of the hypergeometric function 
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we can perform the differentiation with respect to A in Eq. ()B16|) and integrate over (3 to get 
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Similar to the above one can show that the dilogarithms in Eq. ()B18|) are numerically small and 
can be neglected compared to the rest of the expression. The Fourier transforms in the leading 
term in the expansion of running coupling corrections in Eq. ()B11|) was performed in obtaining 
Eq. (j43p (see also the derivation of Eq. ()A12|) ). That result, combined with Eq. (jB18|) . allows 
us to write 
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This expansion again demonstrates the emerging geometric series when higher order fermion 
loops are included in K,^^^. Resumming those series to all orders we obtain the "longitudinal" 
part of Eq. (IU7j) . 

Finally, adding Eq. ()B10|) and Eq. ()B19|) together yields 

ICq{xo, Xi; z) = iC^ixo, Xi, z) + IC^{xo, Xi; z) 

4e-f-2-^ 



a^/C "(a;o,iCi;z)a^/52 In 



R'^{xq,Xi; z) ^j? 



MS, 



X 



4g-5/3-27 
1 - «^/?2 In I I ,0 2 1 + 



4g-5/3-27 

1 I ^31^7^1 + 



(B20) 



which, after resumming the geometric series gives Eq. (j^7j) . as desired. 

C Comparison with dispersive calculation of [37] 

Here we demonstrate explicitly that our all orders expression in Sect. |31 Eqs. ()85|) to 
agree with the results found in [37]. In [37] the running coupling effects were presented in 
terms of transverse and longitudinal Borel functions and with Borel parameter u. With 
our convention for the kernel and the shorthand notation ri = Xq — z, r2 = Xi — z and the 
replacement p q, q q' to match notations in this paper, we quote the expressions of [37] 
as 



q2 q'2 _ ^2 



^a^,y;zB (u, Vi^i^, T^^iy^) = J ^^^^^ 6 

(Cla) 

}C^,y■,.B^{u,n^^^,r,^^^) =et- / e^"^'- e"''''- , . (Clb) 

J l^^TTj q q 

Everything else follows from the definitions for the Borel representation of what is called the 
coupling function i?(riA, r2A) in [37] (expressed in terms of the QCD scale A) which takes the 
form of a sum of transverse and longitudinal contributions 

R{riA, T-sA) = ^ ^ du T{u) j {b'^{u, nfiy^, ra/iMs) + ^^(^' ^i/^HS' ^2/iMs)) 

(C2) 

where for one- loop running as employed in this paper T{u) is to be set to one. The notation for 
the /3-function coefficients is such that /52 = Po/n. All that is left to do to compare ()U1|) to our 
present results is to perform the Borel integral. Up to an overall factor of a^/7r, this amounts 

to replacing the Borel powers f-^) by the corresponding geometric series ^ 



l+£of!iiln 
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for both a = q and a = q' . All that is left to do is to factor out a common denominator 
7 TYT^T 7 ^x^- One finds 



MS / / \ \ MS 



ii^ e ^ / / \ \ ii^ e ^ 

MS / / V \ MS 



Q'q" ^ + ^lnf^)) (l + ^ln(^^^^ 



MS / \ MS 



MS 



and 




(C3b) 



The sum of these contributions is in full agreement with (jH^ as advertised. 
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